Selforganisation and sympathetic cooling of multispecies ensembles in a cavity 
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We predict concurrent selforganisation and cooling of multispecies ensembles of laser-illuminated 
polarisable particles within a high-Q cavity mode. Resonant collective scattering of laser light into 
the cavity creates optical potentials which above a threshold pump power transforms a homogeneous 
particle distribution to a crystalline order for all constituents. Adding extra particles of any mass 
and temperature always lowers the pump power required for selfordering and allows to concurrently 
trap atoms, for which high phase-space densities are readily available, in combination with many 
other kind of atoms, molecules or even polarisable nanoparticles. Collective scattering leads to en- 
ergy exchange between the different species without direct collisional interactions. We analytically 
calculate the threshold condition, energy fluxes and the resulting equilibrium phase-space distribu- 
tions and show that cavity-mediated energy transfer enhances cooling of heavy particles by adding 
light particles forming a cold reservoir. Extensive numerical many-body simulations support the 
results of our kinetic analytic model. 

PACS numbers: 37.30.-K, 37.10.-x, 51.10.+y 



I. INTRODUCTION 

Laser-induced light forces are routinely used to trap 
and manipulate a large class of polarisable particles from 
atoms, molecules up to larger objects as microspheres, 
nanoparticles or even protozoae [1]. Laser cooling, how- 
ever, has so far been limited to a finite class of atomic 
species [2], very few kinds of molecules [3] or individual 
vibration modes of nanomechanical objects [4]. Success- 
ful cooling requires specific setups with well-chosen laser 
frequencies and field configurations so that the number 
of laser-cooled species has only slowly grown in the past 
years [5]. 

In principle selforganisation and cooling by coherent 
light scattering in cavities provides a general alterna- 
tive to trap and cool any kind of polarisable particles, 
which can be injected into an optical resonator [6-8]. In 
practise, however, the required particle phase-space den- 
sities and laser intensities to reach a useful regime so 
far have been achieved only for atomic ensembles [9-11], 
where the theoretical predictions were fully confirmed 
and showed very fast cooling to sub-Doppler tempera- 
tures [12]. However, as the required phase-space densi- 
ties are hard to achieve for molecules [13], we propose 
here to generalise the scheme by introducing ensembles 
of different species and temperatures simultaneously into 
a single optical resonator. We show that under quite 
general conditions all species are simultaneously trapped 
and cooled using only a single laser frequency without 
the need of direct interparticle interaction. As our cen- 
tral claim we predict that the simultaneous presence of 
any other species will always increase the total scattering 
rate for each particle and thus improve the total perfor- 
mance of cooling and trapping for each individual species. 
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As the most interesting case we study a mixture of a pre- 
cooled and dense enough atomic ensemble with a hotter 
and much smaller ensemble of molecules or nanospheres. 
While it would be impossible to reach the selforganisa- 
tion threshold for the latter alone, combined trapping 
and sympathetic cooling can be readily achieved without 
the need of collisions. In fact the different particles could 
be trapped at different locations within the cavity. Over- 
all this provides for a general route to cool new particle 
species and also allows a simple general setup for simulta- 
neous multispecies cooling and trapping without the need 
of a tailored laser for each species. As for a single species, 
multispecies cooling can be significantly improved using 
several cavity modes [7, 14]. Similarly, a combination of 
traditional laser cooling methods for atoms with sympa- 
thetic cavity cooling could be envisaged to enhance the 
performance of the combined system. 

The paper is organised as follows: After introducing 
the basic setup and model equations, we calculate the re- 
quired threshold intensity and joint phase-space density 
to reach selforganisation and enter the superradiant scat- 
tering regime. Then we derive approximate expressions 
for the cooling dynamics and the energy flow between 
the different ensembles as well as the asymptotic equilib- 
rium phase-space distributions. These results are finally 
checked by numerical simulations of the selforganisation 
and cooling dynamics using a particle model. 



II. MODEL 

Let us consider a mixed dilute gas of S species of N s 
(s = l,...,S) polarisable particles each inside a high- 
Q optical resonator. They arc illuminated by a standing 
wave of counterpropagating off-resonant laser beams that 
transversely cross the cavity and are close to resonance 
with a single mode as sketched in fig. 1. Particles within 
the overlap region of pump laser and cavity mode reso- 
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FIG. 1: (Colour online) Setup. A multispecies ensemble of 
particles within a cavity transversely illuminated by a laser 
close to resonance with a single cavity mode. Above thresh- 
old the particles order in a regular pattern which optimise 
scattering into the resonator. 



nantly scatter light from the one into the other. Pump 
laser and cavity light give rise to a dynamical optical 
potential for each species, modifying the particles' distri- 
butions. For simplicity we approximate the pump field 
in the interaction region by a plane standing wave and 
consider particle motion along the cavity axis only. In 
single species calculations this approximation proved suf- 
ficient to explain the essential physics of selforganisation 
and cooling [9, 15, 16]. An almost ideal implementation 
of this model could be realised by confining the multi- 
ple ensembles into one-dimensional tubes created by two 
crossed rctro-rcflcctcd pump laser beams [17]. Extensions 
of the theory presented in this work to 3d-motion and 
spatially dependent mode functions are straight-forward 
and lead to only minor quantitative changes in the ap- 
propriate limits. 

In terms of the effective laser pump amplitude rj s , the 
light shift per photon Uq, s and the cavity field amplitude 
a, the combined optical potential for the particles along 
the cavity axis is given by [15] 



hrj s (a + a*) sin(fcx) + HUo tB \a\ 2 sin 2 (fcir). (1) 



Combining the one-particle position and momentum vari- 
ables z = (x,p) and introducing the one-particle Hamil- 
tonian function 



H s (z, a, a*) 



P 

2m, 



(2) 



the semi-classical model equations [18] in Klimontovich's 
formulation [19] then read (s = 1, . . . , S) 



dfK,s dH s df K , s _ dH s d/ K ,, 
dt dp dx dx dp 



0. 



(3a) 
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(3b) 

Here, /k.s(z, t) is the so-called Klimontovich distribution 

satisfying / K , s (z,0) = N^ 1 £j/=i s ( z ~ z j.) witn { z j s } a 
set of initial phase points, £ denotes white noise modelling 
the fields' quantum fluctuations with (£(£)£*(*')) = 8(t — 
£'), =0. k > designates the cavity decay 

rate and A c = lo p — lu c is the mismatch between the 
pump frequency uj p and the cavity resonance frequency 
u> c . The equations (3) are equivalent to a set of stochastic 
differential equations (SDEs) for the particles' positions 
and momenta Zj 3 (t) and the mode amplitude a(t). 



III. SELFORGANISATION THRESHOLD 

Let us decompose the distributions according to 
/k, s (M) = f s (z,t) + 6f s (z,t) with f s (z,t) = (f K , s (z,t)) 
denoting the average over a statistical ensemble of similar 
initial conditions {zj 3 }, a(0) and the realisations of the 
noise process £. The ensemble-averaged Klimontovich 
distributions, called one-particle distribution functions, 
fulfil 



dfs 
dt 



P df s 
m s dx 



d ($ a ) df s 
dx dp 



/ dS<f> s dSfs 
\ dx dp 



(4) 



These equations are exact but not particularly useful as 
such because they do not form a closed set. However, for 
N s —> oo, statistical correlations become negligible dur- 
ing the initial stage of the time evolution and hence the 
one-particle distributions satisfy Vlasov's equation [20] 



dfs 
dt 



P df s 
m s dx 



d<£ s dfs 
dx dp 



0, 



(5a) 



where $ s = Q s (x, (a), (a*)). In the rest of this work we 
shall for convenience omit the ensemble-average brack- 
ets. These equations together with the average form of 
eq. (3b) for the ensemble-averaged mode amplitude 



(5b) 



represent the essence of the Vlasov kinetic theory of po- 
larisable particles in a resonator describing the initial evo- 
lution purely due to the mean field interaction [21]. 

Note that spatially homogeneous initial distributions 
f s (x,p,0) = fo,s(p) with zero cavity field a(0) = 
are equilibrium states of (5). In any finite ensemble, 
however, density fluctuations cause light scattering and 
the particles experience friction and diffusion. This is 
mathematically described by the correlation term in (4), 
which induces a slow "collisional" evolution of the Vlasov 
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This growth finally ceases and particles and field reach 
a quasistationary selforganised state. As a central result 
of this work let us emphasise here, that the right hand 
side of eq. (6) only depends on cavity parameters and all 
terms in the sum on the left hand side are manifestly pos- 
itive. Hence adding any extra species will always lower 
the power needed to start the sclforganisation process, 
regardless of temperature and polarisability of the addi- 
tional particles. 

For thermal momentum distributions the integrals 
in (6) are unity and the condition gets particularly sim- 
ple. At higher temperatures, where (kmmv s ) 2 ^ k 2 +5 2 , 
the first term in the denominators of eq. (7) can be ne- 
glected and the field amplitude's growth rate is given by 



7 = — hi 
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FIG. 2: (Colour online) Joint selforganisation of two species 
starting from a perturbed uniform state above the instabil- 
ity threshold (6), such that species one is six times critical, 
whereas species two is far below its proper threshold. Fig- 
ure (a) shows the position distributions in the final state, 
(c) and (d) the momentum distributions initially (dashed) 
and after selforganisation (solid) as determined from numer- 
ical simulations of the particle system [22]. The circles show 
the predictions of eq. (9). (b) shows the evolution of the in- 
dividual order parameters 9 S = | J f s sin(kx)dxdp\ and the 
corresponding final values given by the adiabatic theory. Pa- 
rameters: Ni = 10 4 , N 2 = 500, m 2 = 10mi, fc B Ti = 10 i hn, 
k B T 2 = 2.5 x 10 5 hn, r]i = 2.4k, 772 = 27.4k and wr,i = 10~ 2 fi;. 



equilibria (5) towards a new equilibrium. As a cen- 
tral question we now determine the stability of uniform 
equilibria, i.e. whether small perturbations are damped 
or amplified in the course of time. Linearising equa- 
tions (5) around a spatially homogeneous steady state 
and solving the resulting equations according to Lan- 
dau [23], one can find the conditions for a dynamical 
instability under quite general conditions. Obviously, for 
trapping to occur the effective detuning must be nega- 
tive, i.e. S := A c — |^,iV s [/o.s < 0. For convenience 
we rescale the steady-state distributions as fa. s (p) = 
(Lm s ii s ) _1 G s (^-), where v s > is a typical velocity of 
the sth species, L is the cavity length and we assume that 
these distributions decay monotonously with \p\. Then, 
such an equilibrium is unstable if and only if 
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kBT s 
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A glance at this expression shows, that the instability also 
grows at a larger rate the more terms contribute to the 
sum. Hence, both the required power and time needed 
to achieve selforganisation is lowered by combining sev- 
eral species. Obviously, if one can reach the threshold 
with one species, the system certainly still self organises 
if one adds a second species. Let us remark that there 
exists a dynamical instability for positive effective fre- 
quency mismatch S as well, but it is connected to heating 
and does not lead to an ordered distribution. Numerical 
simulations indicate that the quasi-equilibrium state into 
which the system evolves in case of instability is close to 
a BGK solution [24] of (5a) and (5b). A BGK solution is 
a stationary solution of (5) where all one-particle distri- 
butions depend on position and momentum solely via the 
Hamiltonian functions f s (z,t) — F s (H s ) and a(t) = o?o- 
The real- valued functions F s are essentially arbitrary and 
the steady-state mode amplitude «o needs to be self- 
consistently determined from (5b). Let us note that in 
the weak coupling regime, i.e. \J2 S N s Uo, s \ <C the 
single-particle actions (10) are nearly invariant during 
selforganisation for a wide range of parameters and in 
this case it is therefore possible to relate the functions F s 
to the unperturbed uniform and unstable states /o, s (p) 
to obtain the selforganised state f%°(z) as 



(9) 



du > 



h\6\ 



(6) 



where J s = kl s for untrapped and J s = kI s /2 for trapped 
orbits [25]. For an illustration the reader may consult 
fig- 2. 



where P denotes the Cauchy principal value and ks,T s = 
m s v 2 /2. In that case, initial density perturbations am- 
plify and the cavity mode amplitude grows exponentially 
in time at a rate 7 > that solves 



(7- 



u G' s (u)du 

2k B T s 7-co (j/kVsf+U 



j2 jw^g 

s=l 



(7) 



IV. KINETIC THEORY FOR THE COOLING 
AND ENERGY FLUX 

Let us now turn to describe the system evolution be- 
yond eqs. (5) in the weak-coupling limit where 5 « A c . 
To this end we introduce the single- particle action be- 
longing to the instantaneous average potential seen by 
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the sth species, 

I„ = ±i- I ^2m s [H s - $ s (x>)] dx', (10) 

in which Q s (x, a) ~ 2hr) s Re (a) sin(fcx) as a valid ap- 
proximation in this limit. The corresponding angle vari- 
able 9 S can be obtained from the generating function 
S s = ±J X y / 2m s [H s - $ s (x')]dx' as S = §§^. Starting 
from any initial condition, in the long time limit the one- 
particle distribution functions to a good approximation 
become functions of the single-particle Hamiltonians H s 
and thus actions I s alone. Statistical fluctuations slowly 
modify these distributions in such a way, that the system 
evolves towards equilibrium in a sequence of BGK states, 
f s (x,p,t) ~ f a (I a ,t) [26]. Defining 



2tt 



sm(kx)e- in6s 6B a (11) 



as well as 



1 f 2lx BT 



(12) 



as the average variation of the action along a "frozen" 
orbit, we obtain the following nonlinear Fokker-Planck 
equations for the one-particle distributions f s (I s , t) and 
the average mode amplitude a(t) 



2Ll + u ^. = — 
dt s dh dh 



(A s f s + B s ^+Y,C BL [f s J s ,] 



(13) 



a = ——I Yl N ^ [ 9o,s.fs dl s . (14) 

The collision operator, i.e. the r.h.s. of (13), consists of 
two contributions due to the fluctuation and decay of the 
mode amplitude 



A s [f s ] = -m^ s ^ n2|g "- s|2 



' ^ \D(irw a )\* 



(15a) 



n 2 \gn, s \' 



(k 2 + S 2 + n 2 uj 2 ) (15b) 



\D(mu s )\ 2 

and a generalised Balescu-Lenard term [27, 28] 

Cbl[/., fi] = SttWJV,^ ^ "j g "' s| ' x 

x J \g' m /6 (uuj s - mw,') (n^f( - m^/ a )dl/, 

(16) 

where the prime means taking the function at /; = I[ and 
u; s is the nonlinear frequency defined by u s = dH s /dI s . 
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FIG. 3: (Colour online) Simulated (solid) and analytical (cir- 
cles) steady-state momentum distributions for two different 
species m-2 — 40mi below threshold averaged over 250 reali- 
sations. The distribution of the lighter particles (a) is given 
by a g-Gaussian with gi = 1.4. The dashed curves represent 
Gaussians corresponding to (p 2 )/m a = ksT t . As Q2 = 1.01 
the distribution of the heavier particles (b) is indistinguish- 
able from a Gaussian. Parameters: JVi = 300, N2 = 200, 
V^^i = VN2V2 = 800cj r ,i, k = IOOcjr,!, A c = -2.6ojr,i 
and NiUo,i — N2U0.2 — — 0.1wr,i. 



For spatially uniform ensembles the actions reduce to 
I s — > p/k and the expressions for the coefficients (15) 
given in [22] are recovered. In the above expressions, 
D(s) denotes the so-called dispersion relation given by 



D(s) 



(S + nf + S 2 - i4:TThSj2 N sV 2 s J 



n \9n,s\ gj s 

s + inuj s 



(17) 

if Re(s) > 0. Apart from the integral terms in the dis- 
persion relation D and the dependence of the nonlinear 
frequencies lo s and functions g ntS and U s on the com- 
mon mean field amplitude a, the cavity-mediated inter- 
species interaction is described more explicitly by the 
Balescu-Lenard collision operator (16). It has been de- 
rived before by several authors [29, 30] and accounts for 
the interspecies heat flow, which is mediated by the cav- 
ity field. As noted before by these authors Cbl involves 
resonant actions I;,/ s , i.e. orbits such that the condition 
nuJi(Ii) = moj s (I s ) is fulfilled. This term is a possible 
source of sympathetic cooling is such a setup. 
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A. Joint equilibrium states 

The set of possible equilibria of (13) and (14) can 
be divided into two classes: spatially homogeneous with 
vanishing average field and inhomogeneous with nonzero 
photon number. The hrst exist only for effective red 
detuning 26 < — wr. s and are stable only below the 
threshold determined by (6). The corresponding phase- 
space distributions can be explicitly calculated to give 
g-Gaussians 



fs,eq(x,p) ~ exp 



1>~ 



2m a k^T st 



(18) 



where the ^-exponential function is given by exp q (u) 
i 

[1 + (1 — q)u] 1-9 and we have set 



1 



kftT* = h- 



i 2 + 5 2 



4\6\ 



(19) 



Their minimal "thermal" energy thus is given by fee?* = 
hn/2 and is reached if 6 = —k. Interestingly, we have 
a very small but finite minimum kinetic energy in very 
good resonators but of course if K ~ wr ;S quantum ef- 
fects have to be taken into account. As q s — > 1 the 
g-Gaussian becomes an ordinary Gaussian. For an ex- 
ample the reader is referred to fig. 3. It follows that as 
long as the joint equilibrium {f s ,eq(x,p), s = 1 . . . S} is 
Vlasov-stablc according to (6), the individual equilibrium 
distributions are independent of each other because all 
energy exchange currents cease. The cooling rates, how- 
ever, are modified by energy exchange between different 
sorts of particles, described by the generalised Balescu- 
Lenard operator (16) and can considerably shorten the 
time for a given species to reach its steady state. Above 
threshold, stable and strongly trapped equilibria require 
—5 wj iS and exist already if the uniform state is only 
weakly unstable. These are Maxwcll-Boltzmann distri- 
butions 



fs,eq(x,p) ~ exp 

with kinetic temperatures 

T ^ 
KBJkin.s •— 



■"0,3 

\6\ 



(20) 



(21) 



where the trap frequencies are given by Wg s = 
4r] s u>^ s \ Re a\. Again a is determined from (14) using 

9°' s — 2m Jo * n ^ S ^ m ^- Here the equilibria are indeed 
modified by the presence of additional species through 
the selfconsistent cavity field but the mutual interaction 
is not enough to equalise the kinetic temperatures and, 
as in the uniform case, all equilibrium inter-species heat 
fluxes vanish. Figure 4 shows an example of a jointly 
selforganised steady state. Let us note that for any given 
deeply trapped species the equilibrium uncertainty prod- 
uct AxAp = fcB^kin s/^o s 1S bounded from below by h, 
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FIG. 4: (Colour online) Selforganised steady-state momen- 
tum distributions of (a) species one and (b) species two with 
77i2 = 10mi. In (c) we show the time evolution of the ki- 
netic temperatures and the photon number. The dashed lines 
represent the predictions of eq. (20) and eq. (21). The dash- 
dotted line shows the maximally possible photon number. 
The initial rise in the kinetic temperatures originates from the 
fast initial growth of the cavity intensity above threshold. Pa- 
rameters: JVi = 300, N 2 = 200, v^rn = y/Nsip = 600wa,i, 
k = IOOcjr.i, S = —k and E-r i = Sur. i. 



and thus by twice the minimal value for a particle in 
a classical potential. The additional uncertainty may 
therefore be attributed to the quantum fluctuations of 
the mode amplitude. The energy per particle E s can be 
shown to be 



E s = AxApLJo,s > tit*) 



0..s- 



(23) 



AxAp > h, 



(22) 



which is again twice the usual value. The minimal un- 
certainty state, which coincides with the minimal en- 
ergy state, is attained if 2wo,s 3> K for a detuning 
6 = — 2wo, s - These findings remain correct in an en- 
tirely quantum-mechanical treatment because for deeply 
trapped particles and thus approximately harmonic po- 
tentials, the semiclassical equations are exactly equiva- 
lent to the quantum equations. 



B. Sympathetic cooling 

Let us finally examine the energy flow per particle 
Q2->i from species one to species two for two spatially 
homogeneous ensembles a little closer. If we assume that 
species one is already cold, i.e. 2fcBTki n ,i/?iK <C k/wr ; i 



6 



and far from instability, the inter-species heat flow is es- 
timated to be 



TViTOi 
N 2 m 2 



Q 



IM 



fcsTi 



1 — 2 



i 



2 2 40FM 2 



1 



(k 2 + <5 2 ) 2 



toiT 2 



m 2 T\ 



-3/2 



(24) 



where we wrote T s instead of Tkin,s f° r simplicity. Not 
surprisingly it is maximal if 6 = —n. The proportionality 
of the energy flow to the number of cold particles N\ im- 
mediately hints towards a sympathetic cooling scheme, 
in which a cold ensemble is coupled to a smaller num- 
ber of hotter and heavier particles whose cooling rate is 
enhanced due to this exchange current Q 2 ^i- 

From eq. (16) one sees that if at least one of the sys- 
tems is spatially nonuniform (ordered), the inter-species 
heat flows are effectively suppressed due to the loss of res- 
onances and sympathetic cooling gets inefficient. Hence 
in order to get a useful inter-species energy transfer one 
needs to have a large number of cold particles without 
crossing the instability threshold. This behaviour is ex- 
hibited in fig. 5, where we see that the kinetic energy of 
the heavy particles decays slower when selforganisation 
starts and a field is being built up. However, as the field 
and thus the potential grows the spatial confinement of 
the particles is continuously increased. This behaviour 
can be expected to improve using a larger number of 
non-degenerate modes and several pump frequencies. In 
this manner, each mode contributes to the energy trans- 
fer and dissipation without receiving enough scattered 
photons to actually trap particles. 



V. CONCLUSIONS 

We have shown analytically and in simulations that 
the selforganisation threshold condition for an ensemble 
of several different species of particles inside a transver- 
sally pumped standing wave resonator is strictly reduced 
below the value for each of the species separately. Hence 
the joint selforganisation threshold power is the lower 
the more species present. If the threshold can be reached 
with one species, many species can be added and simul- 
taneously trapped and cooled. In the long-time limit the 
achievable temperatures are only limited by the resonator 
linewidth and get close to the quantum ground state for 
deep traps. In the multi-species case the cooling time 
of a given (heavy) species can be reduced due to energy 
exchange with a second already colder (light) species. 
In general one only needs a single intense and narrow 
laser, frequency-stabilised relative to a cavity mode, to 
simultaneously trap and cool a large number of different 
particles within the same volume. The method thus can 
be applied to atomic gas mixtures, atom-molecule mix- 
tures or even microbcads in a dilute atomic gas. Using 
more complex setups involving several laser frequencies 
and modes can be expected to significantly enhance the 
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FIG. 5: (Colour online) Time evolution of the kinetic temper- 
atures for two species, one heavy and the other lightweight. 
Upper plot: Optimal cooling curve of the heavy species alone 
(blue dashed) vs. cooling curve in the presence of a lighter 
species (blue solid). The cooling is more effective due to 
the exchange of energy between the species. Parameters: 
m 2 = 200mi, ATj = 200, N 2 = 200, v^^i = 134wr,i, 
^N 2 "q2 = 134cjr.i, k = 200cjr,i and S — —k. Lower plot: 
the dashed lines represent the temperature evolutions of each 
species in the absence of the other species. The sympathetic 
cooling effect can be observed initially but as the system 
crosses the selforganisation threshold the inter-species heat 
flow ceases. This can be inferred from the kinetic tempera- 
ture curve of the heavier particles (solid blue) becoming par- 
allel to the curve depicting the heavy particles alone (dashed 
blue). Parameters: iVi = 320, N 2 = 500, v^^i = 207wr,i, 
= 258(xjr,i, k = 200cjr,i and S — —k. 



cooling and lead to more complex distributions of the 
particles. The method can be easily generalised to mov- 
ing ensembles in arbitrary cavity geometries, e.g. ring 
resonators. Here the cavity field mediated interaction of 
the ensembles transfers a stopping force applied to one 
ensemble to any other particles. As the general effect 
has been successfully experimentally demonstrated for 
single-species setups [12, 31, 32], we are confident that 
the multispecies generalisation proposed here are well 
within reach of current technology. This might include 
even atomic hydrogen, which could be stopped sympa- 
thetically with a Helium beam. 
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